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ABSTRACT 

The time and frequency spreading of a shallow water acoustic channel can be 
characterized by its scattering function. In order to exploit some of the time 
delay regularities in the acoustic channel, the scattering function is calculated 
using a Pekeris waveguide model. The correlator loss, which is a measure of the 
amount of signal spreading caused by the channel, is calculated as a function 
of range and depth for three well known classes of signals used in active sonar. 
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Scattering Function of Shallow Water Channels 


EXECUTIVE SUMMARY 

Quieter submarines, as a result of improved noise and signature reduction, operating in 
littoral environments pose a challenging ASW detection problem for active sonar. Regard¬ 
less of the quality of the sensors used and the skill of the sonar operators, the environment 
will still impose limitations on the sonar system’s detection ability. One immediate prob¬ 
lem is the lack of detailed knowledge of the environment (as well as the target and other 
operational factors) that can be used to develop suitable new algorithms. However, even 
if this knowledge were available, the amount of processing required would be unlikely to 
be cost effective. 

Accepting that the complexity of the ocean environment can limit the level of acoustical 
information that can be sensed and usefully exploited, it is then reasonable, at least as a 
first step, to model the undersea acoustic channel as a random time-varying linear filter. 
For example, the scattering function, which is a measure of the second order statistics of 
the random time-varying filter, can be used to quantify the channel and hence ultimately 
sonar system performance. 

In this paper the time and frequency spreading of a shallow water acoustic channel 
is characterized by its scattering function. In particular, its time spreading structure is 
calculated by using the Pekeris waveguide model. The Pekeris waveguide with a uniform 
sound speed profile is an attractive model from the viewpoint that only a minimal number 
of parameters are needed, but is sufficiently complex to exhibit many of the salient features 
observed in actual waveguides, at least to first order. By using such a model it is shown that 
signals with large time-bandwidth products are severely distorted by time and frequency 
dispersive channels. Although large time-bandwidth product signals are generally desirable 
for detection, clutter reduction and range resolution, it is clear that when operating in a 
dispersive channel these signals are of limited value when processed using conventional 
techniques such as simple matched filtering. 
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1 Introduction 


Shallow water provides a challenge for designers of active sonar systems due to its fad¬ 
ing multipath structure. The time and frequency spreading caused by the channel can 
adversely distort the transmitted signal and hence result in a degradation in sonar perfor¬ 
mance. In this paper, the limits of coherent pocessing (in the form of replica correlation or 
its matched filter equivalent), due to pulse distortion caused by the underwater medium, 
will be investigated. A simple and reasonably adequate model of an undersea acoustic 
channel is that of a random linear time-varying filter[1, 2, 3, 4]. The random nature is 
due to a lack of detailed knowledge of the acoustic environment (sound speed profile, 
bottom properties, surface motion etc) and the associated system parameters such as the 
geometry of the transmitter, receiver and target. Therefore, one is dealing with a complex 
dynamical system which warrants a stochastic description. For example, the scattering 
function which is a measure of the second order statistics of the random time-varying filter 
will be used to quantify the channel. This, however, is not to suggest that deterministic 
methods, such as wave equation solutions, are not useful. In this paper it is assumed that 
any regularities, as observed in an averaged sense, can be described by the deterministic 
wave equation. In particular, the Pekeris model will be used to account for the time delay 
structure of the scattering function. 

The remainder of this paper is organized as follows: Section 2 provides the background 
for the scattering function description of the channel. Although the theory can be applied 
to the sonar detection and estimation problem in general, this paper will concentrate on 
one way transmission so that complications from target scattering need not be introduced. 
Section 3 gives details of the Pekeris waveguide. Section 4 applies the results of the previous 
sections to (a) a fully analytical case and (b) numerical cases for three classes of well known 
active sonar transmissions. Finally, section 5 concludes the paper. 


2 Scattering function 


Consider an active sonar system in which the transmitter and receiver are in separate 
locations. Let the transmitted signal be given by s ( t ). By modelling the channel as a 
time-varying linear filter, the received signal r (<) can be written as 

r (t) = J h (r, t)s(t — T)d,T + n (t), (1) 

where h(r,t) is the channel impulse response at time t for a given delay r and n (t) is a noise 
term representing sea and thermal noise. For the purposes of this paper the ambient noise 
will not be discussed any further. The random impulse response h (r, t) can considered as a 
member of an ensemble which describes the multipath and Doppler spread of the channel, 
where r, for a fixed t, is a measure of the time delay spread a signal undergoes, and for 
a fixed r, t indicates how rapidly the channel is changing with time, i.e., the amount 
of Doppler spreading. Usually, the channel response is characterized principally by its 
first and second moments, although higher order moments have been considered. The 
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first moment is taken to be zero due to the uncertainty in the phase, leaving the second 
moment, i.e., the covariance function, to carry the bulk of the statistical information about 
the channel. 

In order to proceed further in the simplest way, two common assumptions are made. 
The first is that the channel response is a stationary random process for a given delay 
(range) interval 1 . This implies that the covariance function depends only on the time 
difference between measurements of the channel response. The second assumes that re¬ 
flection or scattering of the channel is uncorrelated as a function of delay (range). These 
two assumptions then enable the covariance function to be written as 

E {h(ji,ti) ,h* (r 2 ,t 2 )} = K ( ti , At) S {t\ - t 2 ) , (2) 

where E{ } denotes ensemble average, At = h — t 2 , and K ( t , At) is a two variable 
function dependent on the properties of the channel. From this two variable function it is 
useful to define the scattering function, a (r, /), such that 

/ OO 

a(rJ)^ At df. (3) 

-OO 


Now assume that the receiver processes the received signal using a matched filter or 
correlator of the form 

y(T,<t>) = j r ( t ) s* ( t — t) e~i 2n<l>t dt, (4) 

where t and (f) are the relative delay and Doppler of the matched filter respectively. By 
using Eqs. (1) and (2) in (4), the mean power output of this matched filter is 

E{\y{T,<f>)\ 2 } = J J a(A,/) |x(A-r,/ -<f>)\ 2 dfd\, (5) 


where 

X (r, <f>) = j s ( t ) s* ( t - t) e-^dt ( 6 ) 

is the ambiguity function of the pulse s ( t ) and may be interpreted as the output of the 
matched filter for a received signal propagating through an ideal distortionless medium. 
Then the scattering function can be thought of as a two dimensional filter modifying the 
ideal channel matched filter response. Note that the ideal channel (i.e., a single path with 
no frequency spreading) is represented by a (r, /) = S (r) 5 (/). 

Consider a time varying multipath channel described by its scattering function. A 
measure of how much the channel spreads the transmitted pulse can quantified by the 
correlator loss, defined as 

CL = 101og 10 (F ), (7) 

where the loss function, F, is given by 


£{|y(o,o)| 2 } 

lx(o,o )| 2 


( 8 ) 


1 This is equivalent to saying that the channel response is uncorrelated as a function of Doppler shift. 
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This is simply comparing the mean response of a time varying multipath channel to the 
ideal channel. It is convenient to normalize the scattering function by the quantity 

I ja(X,f)dfd\ (9) 

and also normalize the waveform, s ( t ), to have unit energy so that x (0,0) = 1. The loss 
function then becomes 

' ' ' 

Since the scattering function can be determined for any range and depth of interest, the 
normalization ensures the loss function is independent of the absolute signal level at any 
given position. 

The main purpose of the loss function is to give an indication of how the energy of the 
transmitted signal is distributed by the channel. For example, values of the loss function 
close to one indicate the existence of a dominant path, while values much less than one 
suggest a number of different possibilites. It could be that the channel supports multiple 
paths or the channel distorts the signal by other means (e.g., frequency spreading) or 
a combination of both. The loss function by itself does not give enough information to 
distinguish between these possibilities. 

It is clear that having some knowledge of the scattering function is essential in under¬ 
standing and perhaps exploiting the effects of the environment on active sonar transmis¬ 
sions. To this end, assume that the scattering function can be written as 

cr (r, <t>)=P (r) S D (<j )), (11) 

where P (r) is the delay power profile and Sd {4>) is the Doppler spread profile. Fur¬ 
thermore, the delay power profile is assumed calculable from the solutions of the wave 
equation subject to appropriate boundary condtions whereas the Doppler profile, without 
any apriori knowledge, is assumed to take the form 



where Bd is the variance of the Doppler spread. Note that the assumed separability of 
the scattering function into the delay power profile and the Doppler spread profile implies 
that the mechanisms responsible for Doppler spreading are statistically independent of 
the range interval under observation. This assumption is also motivated by how one 
would go about measuring the scattering function without sophisticated signal processing. 
To see this, consider a time varying multipath channel. The delay power profile can be 
determined by probing the channel with a short pulse and measuring the average echo 
power as a function of delay. The Doppler profile can be measured by CW (continuous 
wave) testing of the channel to give the frequency spectrum of the average received power. 


3 Pekeris waveguide 


In this section the delay power profile will be calculated by using the Pekeris waveguide[5] 
to model the shallow water channel. The geometry of the channel is shown in Fig. 1. The 
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Pekeris model is chosen because it exhibits many of the characteristics of real waveguides 
without being overly complicated mathematically. The channel is completely specified by 
the depth of water, H, the sound speed profiles of the water and bottom, c(z) and c b (z), 
the densities of the water and bottom, p and p b , and a pressure release at the air-water 
interface. In fact closed form solutions can be found for special cases of interest [6]. The 
pressure time series as a function of range and depth for an arbitrary source waveform[7] 
can be expressed in terms of normal modes 2 , i.e., 


p (<, r, z) = f {-ji r) 52 ^m(zo)tM2)tfo 2) (k m r) S (/) e? 2 * ft df , (13) 

J 771 


where S (/) is the Fourier spectrum of the transmitted waveform, s ( t ), 'tpm( z ) is the mth 
normal mode eigenfunction as a function of depth, Hq ^ is a Hankel function, and k m is 
the horizontal wavenumber which is related to the sound speed profile, c(z), the source 
frequency, /, and vertical wavenumber, j m , by 


*mW+iSi= 



(14) 


In order to probe the channel’s multipath structure, a very broadband source should be 
used. In doing so, the delay power profile can be estimated from Eq.(13) as 


P(r) = 




(15) 


where the Hankel function has been replaced by the exponential in the long range ap¬ 
proximation, to is beginning of the observation time window (or is the time when P (0) is 
maximum) and r is time delay relative to £q- 


Two cases with uniform sound speed profiles will be considered: (i) rigid bottom and 
(ii) fluid bottom. The rigid bottom case can be solved exactly and the delay power profile 
is given by 


m=0 i=l y CI " 


6(r + to- T im ) 


(16) 


where 


cTi m = \Jr 2 + (z — zo + 2 mH) 2 

CT2m = \Jr 2 + (z + zo + 2mH) 2 
CT3m = \/r 2 + (z + zq — 2 (m + 1 )H) 2 

CT im = \Jr 2 + (z-zo- 2(m + l)H) 2 . (17) 


The derivation of this result is given in the Appendix. 

2 Note that representations of solutions to the wave equation other than normal modes can be used. It 
turns out the normal mode method is especially convenient for calculating the sound field as a function of 
range and depth. 
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For the fluid bottom case the delay power profile is calculated by evaluating Eq.(15) nu¬ 
merically because there is no simple analytic solution available. The modal eigenfunctions 
for the water column axe given by 

^m(z) = Am sin (7 m z ), (18) 


where 


A^ = 27 m [7 m H - cos (7 m H) sin (7 m H) - b 2 sin 2 (7 m H) tan (7 m tf)] 1 (19) 

is a normalization constant 3 and b = p/ pb- The corresponding eigenvalues axe obtained 
by solving 


tan (7 m H) = - 


ImH 


( 20 ) 


where 


a = 2ixfH 


!_! 

c 2 cf‘ 


( 21 ) 


Note that the existence of analytic solutions is valuable in that they reveal physical 
insight to the problem and provide a check, in limiting cases, on numerical solutions to 
problems where exact solutions axe either impossible or extremely difficult to obtain. This 
is the role of the rigid bottom case and it is clear that in going from the rigid bottom case 
to the fluid bottom case the computational difficulties have greatly increased. 


4 Applications 


4.1 Analytic example 


Consider a Gaussian LFM signal transmitted through an ideal rigid bottom shallow water 
channel. A complete analytical expression of the loss function is available in this case. Let 
the transmitted signal have a Gaussian envelope of the form 


s 


/ 2 \i 

/t\ 2 ' 


r t 2 ] 

\7rT 2 / eXP 

[-(r)J 

exp 

JkW- 




( 22 ) 


where T is the effective duration of the pulse and W is frequency sweep width (the CW 
case is obtained by setting W = 0). By using Eq.(6), the ambiguity function is[4], 


I = exp 



exp [-7T 2 (4>T - Wt ) 2 ] . 


(23) 


3 The wave solutions of the depth dependent part of the Helmholtz equation form am orthonormal set. 
The normalization constant is obtained by having these solutions satisfy the boundary and orthogonality 
conditions. 
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The loss function of Eq.(10) for the rigid bottom channel ( via Eqs. (11), (12) and (16) ) 
is 

40 * < 24 > 


1 °° ** 1 roo 

F =- J== E E 7- ^2 / lx (Ar im , 0)1 exp 

N^irBjj m=0i=l ( CT *»w) •'~°° 


where Ar, m = — to is the excess delay and the normalization factor is given by 


oo 4 


»=EE 


m=0 i=l ( CT im) 

For the Gaussian pulse, the loss function becomes 

1 ^ ^ 1 


2 • 


oo 4 

EE 


Ny/l + 2ir 2 B 2 D T 2 m=0 ^ ( CT im) 


-2 exp 


Ar im 


7r 2 W 2 Ar 2 m 

T 2 

6Xp 

(1 + 27T 2 B 2 T 2 ) 


(25) 


(26) 


This example explicitly shows how a signal cannot be transmitted with an arbitrarily large 
pulse duration and/or bandwidth in a dispersive channel without incurring significant 
signal spreading. More specifically, F <C 1 if 


T » T c = 


W » B c = 


1 


\Z2ttBjj 

1 


and/or 


it At 


(27) 


where At is the average excess delay, and T c and B c denote the coherence time and 
coherence bandwidth respectively. Conversely, the dispersive effects of the channel can 
be considered unimportant if the sonar waveform’s duration and bandwidth are operated 
within their coherence limits, though this might not be possible in practice. Note that 
these conclusions hold in general though the expressions for T c and B c will differ according 
to the channel and waveforms under consideration. 

Figure 2 shows an example of the correlator loss as a function of range and depth 
assuming a source depth of 10 m and a rigid bottom of depth 100 m. Fig. 3(a) gives an 
example of the Gaussian LFM ambiguity function. It can be seen from Fig. 2 that the 
loss function decreases with range and is largely insensitive to changes in depth. When 
the source and receiver are relatively close to each other the direct path is dominant as the 
other paths have longer path lengths and are weaker in strength (recall that the amplitude 
of each path is inversely proportional to the path length). This means that the loss function 
should be near unity if, for the moment, Doppler spreading is neglected. As the range 
increases, the direct path becomes less dominant relative to the other paths and the spread 
of the delay power profile becomes larger, resulting in a smaller value for the loss function. 
At very long ranges, all the paths are nearly comparable as the path difference between 
paths falls off inversely with range (see Eq. (17)) and hence the loss function decreases 
slowly with range. This is analogous to the well known phenomenon of energy spreading 
in a shallow water channel, i.e., the transition from spherical to cylindrical spreading. 

An order of magnitude estimate of the loss function or the correlator loss can be 
obtained by considering the following artificial, but instructive, example. Assume that 
the received signal from a multipath channel can be represented by a sum of time delayed 
replicas of the transmitted waveform and each path is roughly equivalent. Then the 
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received signal can written as r (t) = 1 /y/N Y,k=i <t>k (*) where N is the number of paths 
and <f>k{t) is the k th delayed replica of the transmitted waveform. The correlator output 
for the received waveform has a series of peaks with amplitude ~ 1/y/N. This gives a 
correlator loss of CL ~ —10 log 10 N. For example, 10 to 100 paths results in a correlator 
loss of —10 to —20 dB, in broad agreement with the results shown in Fig. 2. Also notice 
that the correlator loss for the LFM signal is worse than that of CW pulse. This is not 
unexpected since the larger bandwidth of the LFM signal allows it to resolve a greater 
number of paths. 


4.2 Numerical examples 


In this section the correlator loss as a function of range and depth will be calculated for 
both the rigid bottom and fluid bottom cases using three well known classes of sonar 
signals, as characterized by their ambiguity functions[3, 8]. They are the (i) ridge (e.g., 
rectangular CW pulse), (ii) sheared ridge (e.g., rectangular LFM pulse) and (iii) thumbtack 
(e.g., pseudo-random noise) ambiguity functions given respectively by 


\x{r,4>)\ 2 = 
\x(r,<t>)\ 2 = 
\x(r,(f>)\ 2 = 


+ 


1 1 
s 

4 1 

sin 7 x<j>(T — |r|) 
tt</> (T — |r|) 

", \r\<T 

(28) 

rr-M 

sin7r£ (T — jrj) * 

2 

, M < T 

(29) 

T 

1 

fc 

\T~\t\ 

sin7rrW sin7r</> (T — |r|) ' 

2 


T 

■ktW 7T0(T— |t|) 



TW ( X 


VI 

~ 

VI 

(30) 


where £ = </> — Wt/T and the last equation (representing the thumbtack ambiguity func¬ 
tion) assumes TW > 1. The ambiguity functions axe zero for |r| > T and the last term of 
the thumbtack ambiguity function is zero when \<f>\ >W. The delay-Doppler characteris¬ 
tics of these ambiguity functions are illustrated in Figs. 3(b), 3(c) and 3(d) respectively. 
The following parameter set will be used for the examples: 


Pulse duration, T = 0.5 s 

Frequency sweep width, W = 0 for CW and W = 200 Hz for others 

Observation window length of 2 s 

Centre frequency, f c = 1000 Hz 

Frequency window, 800 - 1200 Hz 

Sampling frequency, f s = 4096 Hz 

Source depth, zq = 10 m 

Channel depth, H = 100 m 
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Sound speed in water, c — 1500 m/s 

Density of water, p = 1000 kg/m 3 

Sound speed of “typical” bottom, q, = 1600 m/s 

Density of “typical” bottom, pb = 1500 kg/m 3 

Sound speed of low reflective bottom, Cb = 1510 m/s 

Density of low reflective bottom, pb = 1010 kg/m 3 

Doppler spread, Bd — 1 Hz 


Note that the Doppler spread 4 is calculated by assuming a range rate spreading of v ~ 2 
knot which gives a Doppler shift of about 1.4 Hz. So Bd — 1 Hz is taken as the nominal 
value. 


Closed form expressions of the loss function for the rigid bottom case using the ridge 
and th umb tack ambiguity functions can be calculated from Eqs.(28), (30) and (24). There 
is no simple form for the sheared ridge ambiguity function but a relatively simple upper 
bound can be obtained. Their expressions are given as follows: For the ridge ambiguity 


function, 


F = 


EE 


NV2^B D T ^0 fri ( CTim) 


f 11 


(31) 


while for sheared ridge 


F < 


EE 


N\/2ttB D T ^ ( CTim) 


/(Ar im ) exp 



WAT im \ 
B d T ) 



(32) 


and for the thumbtack ambiguity function 

1 ^ A 1_ ... f sin 2 (irWArim) Bp \ 

N^B D T^^(cr im ? n T 1 (*WAT im ) 2 W /’ 

where . . , 

[ i _ 111 if 1-7-1 < T* 

f ( T ) = ) 1 T 11 l T l ^ x 

v I 0 otherwise, 


J = \pht erf ^ 


W \ 2 B d 

BdV^J W 


_ w z 

1-e 


(33) 


(34) 


(35) 


k = Tf (Atim) and erf(-) is the error function 5 . The loss functions of the three signal 
types are shown in Figs. 4, 5 and 6 as a function of range and depth. 

4 The Doppler shift is given by <j> = 2f c v/c ~ 0.7 (j^j) (*^t) Hz - 
s The error function is defined as erf(x) = f 0 ex P ( — t 2 ) dt. 
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Unfortunately, no closed form expressions are available for the fluid bottom case. As a 
result the scattering and loss functions are evaluated numerically for each target range and 
depth of interest. An example of the scattering function for a range of 10 km and a depth of 
10 m is shown in Fig. 7. The correlator loss for the three signal types of a “typical” Pekeris 
waveguide is shown in Figs. 8 to 10. Figures 11 to 13 show the case when the waveguide 
has a less reflective bottom. It can be seen that the correlator loss, a quantitative measure 
of signal spreading, decreases in going from a fully reflective bottom to a less reflective one 
for all signal types. Note that the number of paths increases when approaching the rigid 
bottom limit. In such cases the transmitted signal suffers more distortion since more paths 
contribute to the interference. At the same time, the superposition of a large number of 
paths tends to have an averaging or smoothing effect on the correlator loss (see Eq.(10)) 
as a function of range and depth. Furthermore, the signal distortion increases with range 
since the number of significant paths for a given arrival point should also increase with 
range. Finally, since the channel under consideration is predominantly a delay spread one, 
signals with large time-bandwidth products, such as LFM and PRN, will suffer the most 
degradation 6 (see Eq.(27)). 


5 Conclusions 


The scattering function provides a useful statistical characterization of acoustic undersea 
channels in active sonar applications. In particular, the shallow water channel can be 
modelled by the Pekeris waveguide. The Pekeris waveguide with a uniform sound speed 
profile is an attractive model from the point of view of needing only few parameters but 
is otherwise sufficiently complex to exhibit many of the salient features observed in actual 
waveguides, at least to first order. By using such a model it was shown that signals with 
large time-bandwidth products axe severely distorted by dispersive channels. Although 
large time-bandwidth product signals are generally desirable for detection, clutter reduc¬ 
tion and range resolution, it is clear that when operating in a channel with time and 
frequency spreading these signals axe of limited value when processed using conventional 
techniques such as simple matched filtering. 
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6 Notice that for the two CW cases (fluid bottom) in Figs. 8 and 11, there appears a ridge in the 
correlator loss function for depths near the bottom of the channel. This is caused by the consistently 
nonzero value of the modal amplitudes near the bottom of the waveguide (see Figs. 5.3 and 5.15 of 
Ref.[6]). For broadband signals, the sum over a large number of frequency components is expected to wash 
out this effect. A similar argument applies for channels with highly reflective bottoms where the existence 
of a larger number of trapped propagating modes can also cause a cancellation. 
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6 Appendix 


For an isospeed rigid bottom channel the modal eigenvalues and eigenfunctions in Eq.(13) 
are given by 

= ( rn + \') Tl form=:0 ’ 1 > 2 >--- 


7 m 


tpm(z) = Sin ( 7 m z) 

respectively. By using the following results 


^ 2 sin (' y m z 0 ) sin ( 7 m z) H^ ( k m r) = ^ sin ( 7 m z 0 ) sin ( / y m z) ( k m r ) (39) 


7m — 7m—1 

00 


(36) 

(37) 

(38) 


( 2 ) 


m=0 


m=—00 


-jjVc 2 t 2 —r 2 


/ oo 

S(f)^ t df = s(t) 

-OO 


(40) 

(41) 


Eq.(13) becomes 


c 00 f°° 

p(t,r,z) = — s in( 7 m^o)sin( 7 m z) / s(t 

tl * J — OO 


= -31m V c 2 t 2 —r 2 


\/c 2 t 2 — r 2 


-dr 


c f 00 s(t — t) 


s: 


AH J -00 Vc 2 t 2 - r 2 


°° 'I 

jgj7m?l _j_ gj'7m6 _ gi7m?3 _ gjlmU | ^ 


(42) 


where 


c 2 r 2 — r 2 


- z 0 - \/i 
+ z 0 - \/c 2 t 2 - r 2 


£2 = -2 + zo 
£3 = z + z 0 - Vc 2 r 2 — r 2 
£4 = -z - z 0 - \/ c 2 t 2 - r 2 . 

The above expression can be furthered simplified by using 

00 00 

£ ^ T 8 (t — mT) 


8 [f fe)] = 53 


dz 


-1 


£ (a; - Xj), for / (a:,) = 0. 


X=Xi 


The result, after some manipulation, is 


°° f 

P ( t , r, z) = 53 { 

<rr> —_' 


S(t-T a ) S (t Tfe) 


CT„ 


CT h 


}■ 


(43) 

(44) 

(45) 

(46) 


10 



DSTO-TR-0644 


where 

cr a = \jr 2 + (z — zq + 2 mH) 2 

cr b = \Jr 2 + (z + zq + 2 mH) 2 

or in the more commonly used form 


00 4 f—Vu+i 

P(t,r,z) = s{t-T im ) 

m=0 i= 1 ( CTlm > 


(47) 


(48) 


where Ti m are given by Eq.(17). The result of Eq.(16) easily follows from this by using 
s(t) = S (t ). 
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Figure 1: Geometry of the Pekeris waveguide. 
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Figure 2: Correlator loss of (a) a Gaussian CW pulse and (b) its corresponding contour 
plot. Similarly for (c) and (d) but for a Gaussian LFM pulse. The shallow water channel 
is isospeed with a rigid bottom. 
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(a) Gaussian LFM (b) Rectangular CW 



Doppler (Hz) -2 -0.1 Delay (sec) Doppler (Hz) -2 -1 Delay (sec) 
(c) Rectangular LFM (d) PRN 
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Figure 3: Ambiguity functions for (a) Gaussian LFM, (b) rectangular CW, (c) rectangular 
LFM, and (d) PRN (pseudo-random noise). These are shown using T = 0.5 s and W = 40 
Hz (W — 0 for CW) in order to accentuate their features. 
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Figure 4: The correlator loss of a CW pulse (T — 0.5 s) in an isospeed channel with a 
rigid bottom. The contour plot shows four levels: CL = —9.5 dB (solid), —10 dB (dashed), 
—10.5 dB (dotted) and —11 dB (dashdot). 
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Figure 5: The correlator loss of a LFM pulse (T = 0.5 s and W = 200 Hz) in an isospeed 
channel with a rigid bottom. The contour plot shows four levels: CL = —17 dB (solid), 
—18 dB (dashed), —18.5 dB (dotted) and —19 dB (dashdot). 
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PRN pulse - Rigid bottom 
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Figure 6: The correlator loss of a PRN pulse (T = 0.5 s and W = 200 Hz) in an isospeed 
channel with a rigid bottom. The contour plot shows four levels: CL = —18 dB (solid), 
—18.5 dB (dashed), —19 dB (dotted) and —19.5 dB (dashdot). 
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Figure 7: (a) Scattering function and (b) its Delay Power Profile for a target range, r = 10 
km, and depth, z = 50 m. The channel is isospeed with a fluid bottom (c\, = 1600 m/s and 
Pb = 1500 kg/m 3 ). 
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Figure 8: The correlator loss of a CW pulse (T — 0.5 s) in an isospeed channel with a 
fluid bottom (cb = 1600 m/s and pb = 1500 kg/m 3 ). The contour plot shows four levels: 
CL = —3.2 dB (solid), —3.8 dB (dashed), —4.4 dB (dashdot) and —5 dB (dotted). 
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